Signals of the QCD phase transition in core-collapse supernovae 
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We explore the implications of the QCD phase transition during the postbounce evolution of core- 
collapse supernovae. Using the MIT bag model for the description of quark matter and assuming 
small bag constants, we find that the phase transition occurs during the early postbounce accretion 
phase. This stage of the evolution can be simulated with general relativistic three-flavor Boltzmann 
neutrino transport. The phase transition produces a second shock wave that triggers a delayed 
supernova explosion. If such a phase transition happens in a future galactic supernova, its existence 
and properties should become observable as a second peak in the neutrino signal that is accompanied 
by significant changes in the energy of the emitted neutrinos. In contrast to the first neutronization 
burst, this second neutrino burst is dominated by the emission of anti-neutrinos because the electron- 
degeneracy is lifted when the second shock passes through the previously neutronized matter. 

PACS numbers: 12.38.-t 12.38.Mh 21.65.Qr 25.75.Nq 26.30.-k 26.50,+x 26.60.-c 26.60.Kp 95.30.Sf 95.55.Vj 
95.85.Ry 97.60.Bw 



In search of the phase transition from hadronic to de- 
confined matter, heavy ion experiments at RHIC and at 
LHC at CERN explore the QCD phase diagram for large 
temperatures and small baryochemical potentials. For 
these conditions, which were also present in the early uni- 
verse, lattice QCD calculations predict a crossover tran- 
sition between the deconfined chirally symmetric phase 
and the confined phase with broken chiral symmetry. For 
high chemical potentials and low temperatures a first or- 
der chiral phase transition is expected and will be tested 
at the FAIR facility at GSI Darmstadt. 

Due to their large central densities, compact stars can 
also serve as laboratories for nuclear matter beyond sat- 
uration density and may contain quark matter The 
formation of quark matter in compact stars is mainly dis- 
cussed in two scenarios, injprotoneutron stars (PNS) after 
the supernova explosion [2| and in old accreting neutron 
stars (H, 3 ■ For the first case, deleptonization leads to 
the loss of lepton pressure and therefore to an increase 
in the central density so that the phase transition takes 
place. Possible observables are the emission of gravita- 
tional waves [H, 0] due to the contraction of the neutron 
star or delayed 7-ray bursts [B[ . 

In this article we want to follow a third and less dis- 
cussed case. The phase transition from hadronic to quark 
matter can already occur in the early postbounce phase 
of a core-collapse supernova @, 0, II, 0- This re- 
quires a phase transition onset close to saturation den- 
sity, which can be realized for high temperatures and low 
proton fractions. For such a scenario Ref. [|| found the 
formation of a second shock as a direct consequence of the 
phase transition. However, the lack of neutrino transport 
in their model allowed them to investigate the dynamics 



only for a few ms after bounce. Very recently, a quark 
matter phase transition has been considered with Boltz- 
mann neutrino transport for a 100 M Q progenitor 
The appearance of quark matter shortened the time until 
black hole formation due to the softening of the equation 
of state (EoS), but did not lead to the launch of a second 
shock. 

In our core-collapse simulations of low- and 
intermediate-mass progenitor stars, we confirm the 
formation of a second shock caused by the phase tran- 
sition to quark matter. We even find that the second 
shock triggers a delayed supernova explosion during 
the postbounce accretion phase. This represents an 
interesting addition to currently discussed supernova 
explosion mechanisms, such as the neutrino-driven 
the magneto- rotational [13 . fll ] or the acoustic 



mechanisms [15l ] . A clear imprint of the phase transition 
to quark matter at the launch of the explosion could 
be expected in the neutrino signal, the emission of 
gravitational waves and the nucleosynthesis yields. 
Given the improvements in the sensitivity of neutrino 
detectors since SN1987A, it seems that the neutrino 
signal from a Galactic supernova is the most likely 
observable which could provide a first indication of the 
early phase transition to quark matter. 

Lattice QCD is not yet able to make predictions for 
large chemical potentials relevant for neutron star calcu- 
lations. Consequently, the quark matter EoS is currently 
computed using phenomenological descriptions such as 
the MIT bag or the NJL models. As we aim to study 
the basic effects from quark matter phase transitions on 
core-collapse supernovae, we chose the very simple but 
widely applied MIT bag model. In modeling the phase 
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transition to quark matter there is a main physical un- 
certainty: the critical density n cr i t for the onset of the 
mixed phase. In our model, n crit is determined by the 
bag constant B and the strange quark mass. At present, 
the bag constant is not a fixed parameter and is usu- 
ally chosen between S 1 / 4 = 145 - 200 MeV 0. We 
choose the bag constant such that we obtain an early 
onset for the phase transition and a maximum mass of 
more than 1.44 Mq, without enabling absolutely stable 
strange quark matter. Within this narrow range we se- 
lect S 1 / 4 = 162 McV (eosl) and 165 McV (eos2), and 
a strange quark mass of 100 MeV as indicated by the 
Particle Data Group - The two choices of the bag 
constants lead to critical densities of ~ 0.12 fm -3 and 
~ 0.16 fm -3 , respectively (for T = and proton fraction 
Y p = 0.3). For the hadronic EoS we use the table of Shcn 
et al. [10| . The phase transition to quark matter is mod- 
eled by a Gibbs construction as discussed in Refs. [{J Hlj 
and is therefore of second order according to the Ehren- 
fest classification but a phase transition of the first kind 
after Landau's definition The phase diagram using 
eosl for different proton fractions Y p is shown in Fig. [T] 
We remark that larger values for the bag constant result 
in higher critical densities. For the sake of simplicity, we 
have neglected finite size effects and Coulomb interac- 
tions within the mixed phase. Such contributions would 
decrease the density range of the mixed phase making the 
phase transition more similar to a Maxwell construction 
[201 ] . For comparison, we have also computed the EoS by 
using the Maxwell construction and find that the results 
of our core-collapse simulations arc qualitatively similar. 

We note that the small obtained values for the crit- 
ical density close to saturation density are not in con- 
tradiction with heavy ion data. In contrast to heavy- 
ion collisions high-density supernova matter is isospin- 
asymmetric with a proton fraction Y p ~ 0.3. Further- 
more SN timescales are in the range of ms and there- 
fore long enough to establish equilibrium with respect 
to weak interactions that change the strangeness on 
the timescales of micro-seconds or less. The additional 
strange quark degree of freedom and the large asymme- 
try energy allow one to obtain small values for n cr n and 
lead to an early appearance of quark matter (see Fig.]]}. 

Our choice of parameters is also compatible with the 
still most precise neutron star mass measurement of 1.44 
M Q for the Hulse- Taylor pulsar [2l|. With eosl and 
eos2, we obtain values for the maximum gravitational 
mass of 1.56 and 1.50 M© respectively. Note that higher 
neutron star masses can be achieved with more sophis- 
ticated models of quark matter [23] - For _B 1//4 =162 and 
165 McV, almost the entire star is composed of quark 
matter, surrounded by a mixed quark-hadronic phase, 
which is enclosed by a thin pure hadronic crust. 

For the accurate prediction of the three-flavor neu- 
trino signal, general relativistic effects may be impor- 
tant. Hence, we choose for our investigation a well-tested 
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FIG. 1: The onset (thick lines) and the end (thin lines) of 
the mixed phase from the QCD phase transition from nuclear 
matter to quark matter using eosl. 



general relativistic description of the neutrino radiation 
hydrodynamics in spherical symmetry, that is based on 



25| . Quark matter 



Boltzmann neutrino transport 
appears only in optically thick regimes where neutrinos 
are in thermal and chemical equilibrium with matter. For 
this first proof-of- principle study, we use hadronic weak 
interaction rates [261 ] in the quark phase and derive the 
hadronic chemical potentials from the quark chemical po- 
tentials such that weak equilibrium for neutrinos in quark 
matter is obtained. Our simulations are launched from a 
10 and a 15 M Q progenitor model from Ref. 
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The standard core-collapse scenario leads to core 
bounce at nuclear saturation density and the formation 
of a shock. This expanding shock looses energy due to the 
dissociation of nuclei and the emission of the i/ e -burst at 
~ 10 ms after bounce (see Fig. [3] (a)) and therefore turns 
into a standing accretion shock (SAS). The SAS could 
be revived by neutrino heating [l2j. However, explosions 
in spherically symmetric models with accurate neutrino 
transport have only been obtained for a 8 M Q ONcMg 
progenitor star [28|. The collapse of more massive pro- 
genitors leads to an extended postbounce phase, during 
which the central PNS contracts due to mass accretion. 

In our models that allow a transition to quark matter, 
the onset of the mixed phase (rag ~ 0.1 fm -3 , T ~ 10 
McV, Y e ~ 0.3) is already achieved at core bounce. The 
initially reached quark matter fraction at the center of the 
PNS remains small during the first 50 ms after bounce. 
In the subsequent compression, the quark matter fraction 
rises again and an increasing central region of the PNS 
enters the mixed phase. The reduced adiabatic index 
causes the PNS to collapse. 

PNS collapse: The central density rises up to 4 - 5 
times nuclear saturation density where the collapse halts 
due to the stiffening of the EoS in the pure quark phase. 
A large fraction of the PNS is composed of quarks, en- 
closed by a mixed hadronic-quark phase, which is sur- 
rounded by the infalling hadronic envelope (see Table [I]). 
The mixed phase region shrinks gradually during the 
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FIG. 2: Neutrino luminosities (a) and (b) and rms-energies (c) 
measured at 500 km distance for a 10 Mq progenitor model. 
The results of the quark EoS eosl (thin lines) are compared to 
the results of the pure hadronic EoS [HI (thick lines). A sec- 
ond neutrino burst is clearly visible at ~ 260 ms after bounce. 
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"moment of black hole formation 

'black hole formation before positive explosion energy is achieved 

TABLE I: Selected properties of the different models. t p t 
is the time after bounce at the first appearance of the pure 
quark phase. The different mass distributions (baryon mass) 
of the PNS, the quark mass Mq, the mass of the mixed phase 
M m i lel ! and the total PNS mass Mpns are taken at a later 
stage when clearly positive explosion energies E exp i are ob- 
tained. 



PNS collapse as more and more matter converts from the 
mixed into the pure quark phase. On this short time scale 
of ~ 1 ms, the SAS remains almost unaffected by this dy- 
namical evolution inside the PNS (sec Fig. |3J). However, 
the change in the chemical potentials and the increasing 



density during the phase transition establish weak equi- 
librium at a lower electron fraction of Ye < 0.1, while the 
mean energy of the trapped v e increases above 200 MeV. 

Shock formation and early shock propagation: A sub- 
sonic accretion front forms at the interface between the 
hydrostatic pure quark phase and the infalling mixed 
phase (thick dashed line Fig. [3]). The accretion front 
propagates through the mixed phase, meets the super- 
sonically infalling hadrons at the sonic point and turns 
into an accretion shock (thick dash-dotted line). The 
internal energy of the infalling hadronic matter, that is 
accumulated on the compactified core of the PNS, in- 
creases significantly due to the compression in the strong 
gravitational field. This leads to a rapid conversion of 
hadronic matter into the mixed phase. As the accreted 
layers become less dense, the second accretion shock fi- 
nally detaches from the mixed phase boundary and prop- 
agates into the pure hadronic phase. This phase was 
dclcptonized by the continued emission of electron neutri- 
nos after the first neutronization burst. In this electron- 
degenerate environment, weak equilibrium is achieved at 
an electron fraction ~ 0.1. When the second shock runs 
across this matter, the electron-degeneracy is lifted by 
shock-heating and the weak equilibrium is restored at 
higher values of the electron fraction (Y e > 0.2). 

Explosion: As the second shock propagates through 
the outer layers of the PNS, the ram pressure ahead of the 
shock decreases rapidly due to the large decline of density 
and velocity. The pressure behind the shock is supported 
by the accumulation of the infalling shock-heated matter 
and the increasing electron pressure due to the larger Y e 
in the adjusting weak equilibrium. Hence, this accretion 
shock accelerates and turns into a dynamic shock with 
positive matter velocities (thin solid line in Fig. [3]). Up to 
this point, neutrino transport plays a negligible role since 
neutrinos are trapped. This changes when the second 
shock reaches the neutrino spheres. A second neutrino 
burst of all neutrino flavors is released (see Fig. domi- 
nated by v e due to the above-mentioned increase in Y e . In 
addition, the more compact PNS releases (///t) - neutri- 
nos with significantly higher mean energies as illustrated 
in Fig. [2] (c). As soon as the expanding second shock 
merges with the outer SAS, the scenario resembles the 
situation of a neutrino-driven explosion mechanism (thin 
dashed line in Fig. [3J, except for the large matter out- 
flow with velocities ~ 10 5 km/s. Behind the expanding 
matter, a region with matter inflow develops due to neu- 
trino cooling (thin dash-dotted line). The corresponding 
accretion luminosity from the later fallback explains the 
transient increase of the electron neutrino flavor lumi- 
nosities in Fig. [2] (a) ~ 340 ms after bounce. This matter 
inflow in the cooling region becomes supersonic and de- 
velops another standing accretion shock at the surface of 
the PNS at a radius of ~ 50 km. The neutrinos emitted 
from this cooling region are partly absorbed behind the 
expanding shock at a few 100 km radius and support the 
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FIG. 3: Velocity profiles at different times during the post- 
bounce evolution of a 10 Mq progenitor model based on eosl, 
illustrating the development of the explosion through different 
stages; thick solid line at 240.5 ms, dashed at 255.2 ms, dash- 
dotted at 255.5 ms and thin solid line at 255.6 ms, dashed at 
256.4 ms, dash-dotted at 258.3 ms. 



explosion additionally. We obtain explosion energies of 
several 10 50 erg (see Tabic H]). The neutrino luminosities 
decrease after the onset of the explosion (see Fig. [2] at ~ 
350 ms after bounce). 

In general, the models with eosl and eos2 evolve in a 
qualitatively similar manner. However, the models with 
the larger bag constant show a longer PNS accretion time 
before the onset of the phase transition due to the larger 
critical density. This results in a more massive PNS with 
a deeper gravitational potential, so that the second shock 
develops larger explosion energies (see Table H|. In com- 
parison to the simulations using eosl, the second neutrino 
burst appears several 100 ms later and is found to have 
a larger peak-luminosity due to higher temperatures of 
the shocked material. The more massive progenitor stars 
give an earlier onset of the phase transition and result 
in a more massive PNS. A special case is the dynamical 
evolution of the PNS of the 15 M Q progenitor model us- 
ing eos2. Almost simultaneously with the formation of 
the second shock, the more compact quark core collapses 
to a black hole. Still, this does not necessarily exclude 
an explosion. Unfortunately, our co-moving coordinate 
choice does not allow us to follow the dynamical evolu- 
tion beyond the formation of the apparent horizon [25| . 

The main result of this investigation is a strong sig- 
nature of the formation of quark matter in the early 
postbounce phase of core collapse supernovae. A sec- 
ond shock forms inside the PNS, that affects significantly 
the properties of the emitted neutrinos. For a Galactic 
core-collapse supernova, a second neutrino burst should 
be resolvable by the present neutrino detectors. Unfor- 
tunately, the time sequence of the neutrino events from 
SN1987a [3(| was statistically not significant. Anyway, 
the goal of this first study is to predict the general ef- 
fects of the early phase transition to quark matter during 
the postbounce phase. To optimally reproduce the ob- 



servations from SN1987A further analysis and improve- 
ments of the EoS would be required. The magnitude and 
the time delay of the second neutrino burst provide cor- 
related information about the critical density, the EoS 
in different phases and the progenitor model. For low 
and intermediate mass progenitor models, the energy of 
the second shock becomes sufficient to drive an explosion 
even in spherical symmetry. The explosion is powered 
by the accretion of matter into the deeper gravitational 
potential of the more compact PNS, if quark matter is 
present in the core. The ejecta contain neutron-rich ma- 
terial that expands on a fast timescale and should be 
investigated as a possible site for the r-process. With re- 
spect to the remnant, the narrow range of PNS masses 
found in Table [I] may provide an explanation for the clus- 
tering of the observed neutron star masses (gravitational) 



around 1.4 M Q (see e.g. |21(). The discussed direct black 



hole formation at the phase transition could be investi- 
gated further in light of the observed connection between 
supernovae and 7-ray bursts [3ll ]. For other progenitor 
models, a delayed collapse to a black hole may occur after 
the explosion during the PNS cooling phase. 

The presented analysis should be complemented by 
multi-dimensional simulations, to explore the impact of 
known fluid instabilities that can not be treated in spher- 
ical symmetry. Another interesting scenario would be a 
weak neutrino driven explosion, followed by a fallback- 
induced QCD phase transition. Since the QCD phase 
diagram shows a large variety of color-superconducting 
phases 32J, l33l . [34j, a more sophisticated quark matter 
EoS should be adopted. This could lead to a second 
phase transition within the quark core of the PNS and 
would be an interesting extension of the present study. 

This work has been supported by the Swiss National 
Science Foundation under the grant numbers PP002- 
106627/1 and PP200020-105328/1, the Helmholtz Re- 
search School for Quark Matter Studies, the Italian Na- 
tional Institute for Nuclear Physics, the Frankfurt Insti- 
tute for Advanced Studies and by the ESF CompStar 
program. A.M. is supported at the Oak Ridge National 
Laboratory, which is managed by UT-Battellc, LLC for 
the U.S. Department of Energy under contract DE-AC05- 
00OR22725. We would like to thank for stimulating dis- 
cussions with D. Blaschke, A. Drago, G. Martmez-Pinedo 
and S. C. Whitchouse. 



[1] F. Weber, Progress in Particle and Nuclear Physics 54, 
193 (2005). 

[2] J. A. Pons, A. W. Steiner, M. Prakash and J. M. Lat- 

timer, Phys. Rev. Lett. 86, 5223 (2001). 
[3] L.-M. Lin, K. S. Cheng, M. C. Chu and W. M. Suen, 

Astrophys. J. 639, 382 (2006). 
[4] E. B. Abdikamalov, H. Dimmelmeier, L. Rezzolla and 

J. C. Miller (2008), astro-ph/0806.1700. 



5 



[5] Z. Berezhiani, I. Bombaci, A. Drago, F. Frontera and 

A. Lavagno, Astrophys. J. 586, 1250 (2003). [20] 
M. Takahara and K. Sato, Astrophys. J. 335, 301 (1988). 
M. Takahara and K. Sato, Progress of Theoretical [21 
Physics 80, 861 (1988). [22 
N. A. Gentile, M. B. Aufderheide, G. J. Mathews, F. D. [23 
Swesty, and G. M. Fuller, Astrophys. J. 414, 701 (1993). 
A. Drago and U. Tambini, J. Phys. G25, 971 (1999). [24 
N. Yasutake, K. Kotake, M.-a. Hashimoto, and S. Ya- 
mada, Phys. Rev. D75, 084012 (2007). [25 
K. Nakazato, K. Sumiyoshi, and S. Yamada, Phys. Rev. 
D77, 103006 (2008). 

H. A. Bethe and J. R. Wilson, Astrophys. J. 295, 14 [26 
(1985). [27 
J. M. LeBlanc and J. R. Wilson, Astrophys. J. 161, 541 

(1970) . [28 

G. S. Bisnovatyi-Kogan, Soviet Astronomy 14, 652 

(1971) . [29 
A. Burrows, E. Livne, L. Dessart, C. D. Ott and J. Mur- 
phy, New Astronomy Review 50, 487 (2006). [30 
K. Schertler, C. Greiner, J. Schaffher-Bielich, and 
M. H. Thoma, Nucl. Phys. A677, 463 (2000), astro- 
ph/0001467. [31 
S. Eidelman et al. (Particle Data Group), Phys. Lett. [32 
B592, 1 (2004). 

H. Shen, H. Toki, K. Oyamatsu and K. Sumiyoshi, Nucl. [33 
Phys. A637, 435 (1998). 

L. D. Landau and E. M. Lifshitz, Statistical physics. [34 
Pt.l (Course of theoretical physics - Pergamon Interna- 
tional Library of Science, Technology, Engineering and 
Social Studies, Oxford: Pergamon Press, and Reading: 



Addison- Wesley, — cl969, 2nd rev. - enlarg.ed., 1969). 
D. N. Voskresensky, M. Yasuhira and T. Tatsumi, Nucl. 
Phys. A723, 291 (2003). 

J. M. Lattimer and M. Prakash, Science 304, 536 (2004). 
M. Alford et al., Nature 445, E7 (2007). 
A. Mezzacappa and S. W. Bruenn, Astrophys. J. 410, 
740 (1993). 

A. Mezzacappa and O. E. B. Messer, Journal of Compu- 
tational and Applied Mathematics 109, 281 (1999). 
M. Liebendorfer, O. E. B. Messer, A. Mezzacappa, S. W. 
Bruenn, C. Y. Cardall, and F.-K. Thielemann, Astro- 
phys. J. Supplement 150, 263 (2004). 
S. W. Bruenn, Astrophys. J. Supplement 58, 771 (1985). 
S. E. Woosley, A. Heger and T. A. Weaver, Reviews of 
Modern Physics 74, 1015 (2002). 

F. S. Kitaura, H.-T. Janka, and W. Hillebrandt, A&A 
450, 345 (2006). 

A. Steiner, M. Prakash and J. M. Lattimer, Phys. Lett. 
B486, 239 (2000). 

K. S. Hirata, T. Kajita, M. Koshiba, M. Nakahata, 

Y. Oyama, N. Sato, A. Suzuki, M. Takita, Y. Totsuka, 

T. Kifune, et al., Phys. Rev. D 38, 448 (1988). 

T. Piran, Reviews of Modern Physics 76, 1143 (2005). 

S. B. Riister, V. Werth, M. Buballa, I. A. Shovkovy, and 

D. H. Rischke, Phys. Rev. D 72, 034004 (2005). 

D. Blaschke, S. Fredriksson, H. Grigorian, A. M. Ozta§, 

and F. Sandin, Phys. Rev. D 72, 065020 (2005). 

F. Sandin and D. Blaschke, Phys. Rev. D75, 125013 

(2007), astro-ph/0701772. 



